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The most frequently applied algorithm for QR decomposition uses the 
Householder transformation u = Hv , where the Householder matrix H is a 
symmetric and orthogonal matrix of the form: 

H - I - 2xx T 

with the identity matrix J and any normalized vector x with 
z\\l — x?z = 1. 



Householder transformations zero the m — 1 elements of a column vector v 
below the first element: 
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with / = ; 



v /2c(c - v-i) 



fulfils x T z = 1 and that with H = I — xx T one obtains the vector 
[c 0 •■• 0] T . 

To perform the decomposition of the m % h matrix 4 _ qr ( W j tn m > n ) 
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we construct in this way an m x m matrix fit 1 ) to zero the m — 1 
elements of the first column . An ro — 1 x m — 1 matrix will zero the 
m*- 2 elements of the second column. With Qi 2 ) we produce the m x m 
matrix 
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etc 



After 7i (n — 1 for m — n) such orthogonal transforms we obtain: 

R = H (7 ^ k-H^H^A 
R is upper triangular and the orthogonal matrix Q becomes: 

In practice the are never explicitely computed. 
References 
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See Also: Gram-Schmidt orthogonalization 



Other names: Householder reflection, Householder matrix 
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Keywords: matrix orthogonaliz'ation 

Cross-references: upper triangular, orthogonal, column, matrix, decomposition, column 
vector, vector, identity matrix, orthogonal matrix, symmetric, QR decomposition 
There are 5 references to this object. 

This rs version 3 of Householder transformation, born on 2002-01-04, modified 2002-03-08. 
Object id is 1210, canonical name is HouseholderTransformation. 
Accessed 4858 times total. 

Classification: 

AMS MSC: 65F25 (Numerical analysis :: Numerical linear algebra :: Orthogonalization) 



I 5A57 (Linear and multilinear algebra; matrix theory :: Other types of matrices ) 
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Cholesky Factorization 



Cholesky factorization factors an N x N 9 symmetric positive-definite matrix A into the product of a lower triangular 
matrix L and its transpose, i.e., A = LL T (or A = 1/ U, where U is upper triangular). It is assumed that the lower 
triangular portion of A is stored in the lower triangle of a two-dimensional array and that the computed elements of L 

overwrite the given elements of A. At the fr-th step, we partition the n x n matrices A^ k \ L^ k \ and Z/** , and write the 
system as 




where the block ^li is n b * ^21 is ( ra _ "fc) x "a, and ^22 is ( n " n b) * ( n " n b) £11 and ^22 are lower 
triangular. 

The block-partitioned form of Cholesky factorization may be inferred inductively as follows. If we assume that ^n, the 
lower triangular Cholesky factor of -^n, is known, we can rearrange the block equations, 

^21 «- ^2l(£fl) -1 i 

A& <- A& — £21 £31 = £22^22- 

A snapshot of the block Cholesky factorization algorithm in Figure 5 shows how the column panel (^11 and ^21) is 
computed and how the trailing submatrix ^22 is updated^ The factorization can be done by recursively applying the 
steps outlined above to the ( n n *>) x (n — n b ) ma trix ^22. 




Figure 5: A snapshot of block Cholesky factorization. 



In the right-looking version of the LAPACK routine, the computation of the above steps involves the following 
operations: 
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-An ^n^n 



2. dtrsm: Compute the column panel ^21, 



3. dsyrk: Update the rest of the matrix, 




The parallel implementation of the corresponding ScaLAPACK routine, pdpotrf, proceeds as follows: 

1 . pdpotf2 : The process P i, which has the n b x n b diagonal block -Ail, performs the Cholesky factorization of -An. 

o Pi performs -An ~* ^li^fi, and sets a flag if -An is not positive definite. 

o P% broadcasts the flag to all other processes so that the computation can be stopped if -An is not positive 
definite. 

2. pdtrsm; £n is broadcast columnwise by Pi down all rows in the current column of processes, which computes 
the column of blocks of ^21. 

3. pdsyrk: the column of blocks ^21 is broadcast rowwise across all columns of processes and then transposed. 
Now, processes have their own portions of £21 and ^21. They update their local portions of the matrix -A22. 
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